Thermal anomalies associated with shallow gas hydrates in the K-2 mud volcano, Lake Baikal

Thermal measurements and hydrate mapping in the vicinity of the K-2 mud volcano in Lake Baikal have revealed a particular type of association of thermal anomalies (29–121?mW?m^–2) near hydrate-forming layers. Detailed coring within K-2 showed that hydrates are restricted to two distinct zones at sub-bottom depths exceeding 70–300?cm. Temperature data from stations with hydrate recovery and degassing features all display low thermal gradients. Otherwise, the thermal gradients within the mud volcano are generally increased. These findings imply a more complicated thermal regime than often assumed for mud volcanoes, with important roles for both fluids and hydrates. The coexistence of neighbouring low and high thermal anomalies is interpreted to result from discharging and recharging fluid activity, rather than hydrate thermodynamics. It is suggested that hydrates play a key role in controlling the fluid circulation pattern at an early stage. At a later stage, the inflow of undersaturated lake water would favour the dissolution of structure I hydrates and the formation of structure II hydrates, the latter having been observed on top of structure I hydrates in the K-2 mud volcano.     

Quasi-stationary striations in basin-scale oceanic circulation: vorticity balance from observations and eddy-resolving model

Recently, prominent jet-like features of the ocean circulation, called striations, with a meridional scale of O(300–500 ;km) and extending for thousands of kilometers in length, have been detected in satellite and in situ observations and in high-resolution numerical models. In this paper, we study quasi-stationary striations, which are best seen in the multi-year time-averaged velocity fields. Using 1993–2002 mean dynamic ocean topography, satellite altimeter observations, and output of the Ocean General Circulation Model for the Earth Simulator, the dynamics of the quasi-stationary striations in the eastern parts of the subtropical North and South Pacific are examined by assessing individual terms in the time-averaged equations of relative and potential vorticity. While non-linear effects are found to be essential in the dynamics of the striations, rejecting some linear hypotheses forwarded in the earlier studies, the relevance of the Rhines mechanism is not confirmed. Eddy flux does not act as a relative vorticity source for the striations. Using the potential vorticity (PV) diagnostics, we show that the time-mean PV is not conserved along the time-mean streamlines, and on the scale of the striations these changes in PV are largely induced by the eddy flux of layer thickness. The fact that eddy fluxes contribute to the striations’ time-mean PV budget suggests that the striations are not a kinematical artifact of time-averaging of westward-propagating eddies.     

Experimental and petrological constraints on local-scale interaction of biotite-amphibole gneiss with H2O-CO2-（K, Na）Cl fluids at middle-crustal conditions： Example from the Limpopo Complex, South Africa

Reaction textures and fluid inclusions in the～2.0 Ga pyroxene-bearing dehydration zones within the Sand River biotite-hornblende orthogneisses（Central Zone of the Limpopo Complex） suggest that the formation of these zones is a result of close interplay between dehydration process along ductile shear zones triggered by H₂O-CO₂-salt fluids at 750—800℃and 5.5—6.2 kbar.partial melting,and later exsolution of residual brine and H₂O-CO₂ fluids during melt crystallization at 650—700℃.These processes caused local variations of water and alkali activity in the fluids,resulting in various mineral assemblages within the dehydration zone.The petrological observations are substantiated by experiments on the interaction of the Sand River gneiss with the H₂O-CO-2-（K,Na）Cl fluids at 750 and 800℃and 5.5 kbar.It follows that the interaction of biotite-amphibole gneiss with H₂O-CO₂-（K.Na）Cl fluids is accompanied by partial melting at 750—800℃.Orthopyroxene-bearing assemblages are characteristic for temperature 800℃and are stable in equilibrium with fluids with low salt concentrations,while salt-rich fluids produce clinopyroxene-bearing assemblages.These observations arc in good agreement with the petrological data on the dehydration zones within the Sand River orthoeneisses.     

Dynamics of hourly sea level at Hillarys Boat Harbour,Western Australia: a chaos theory perspective

Water level forecasting using recorded time series can provide a local modelling capability to facilitate local proactive management practices. To this end, hourly sea water level time series are investigated. The records collected at the Hillarys Boat Harbour, Western Australia, are investigated over the period of 2000 and 2002. Two modelling techniques are employed: low-dimensional dynamic model, known as the deterministic chaos theory, and genetic programming, GP. The phase space, which describes the evolution of the behaviour of a nonlinear system in time, was reconstructed using the delay-embedding theorem suggested by Takens. The presence of chaotic signals in the data was identified by the phase space reconstruction and correlation dimension methods, and also the predictability into the future was calculated by the largest Lyapunov exponent to be 437 h or 18 days into the future. The intercomparison of results of the local prediction and GP models shows that for this site-specific dataset, the local prediction model has a slight edge over GP. However, rather than recommending one technique over another, the paper promotes a pluralistic modelling culture, whereby different techniques should be tested to gain a specific insight from each of the models. This would enable a consensus to be drawn from a set of results rather than ignoring the individual insights provided by each model.     

A fragmentation criterion for highly viscous bubbly magmas estimated from shock tube experiments

When a highly viscous bubbly magma is sufficiently decompressed, layer-by-layer fracturing propagates through the magma at a certain speed (fragmentation speed). On the basis of a recent shock tube theory by Koyaguchi and Mitani [Koyaguchi, T., Mitani, N. K., 2005. A theoretical model for fragmentation of viscous bubbly magmas in shock tubes. Journal of Geophysical Research 110 (B10), B10202. doi:10.1029/2004JB003513.], gas overpressures at the fragmentation surface are estimated from experimental data on fragmentation speed in shock tube experiments for natural volcanic rocks with various porosities. The results show that gas overpressure at the fragmentation surface increases as initial sample pressure increases and sample porosity decreases. We propose a new fragmentation criterion to explain the relationship between the gas overpressure at the fragmentation surface, the initial pressure and the porosity. Our criterion is based on the idea that total fragmentation of highly viscous bubbly magmas occurs when the tensile stress at the midpoint between bubbles exceeds a critical value. We obtain satisfactory agreement between our simulation and experiment when we assume that the critical value is inversely proportional to the square root of bubble wall thickness. This fragmentation criterion suggests that long micro-cracks or equivalent flaws (e.g., irregular-shaped bubbles) that reach the midpoints between bubbles are a dominant factor to determine the bulk strength of the bubbly magma.     

The Polar Cap (PC) indices: Relations to solar wind parameters and global magnetic activity

The polar geomagnetic activity resulting from solar wind–magnetosphere interactions can be characterized the Polar Cap (PC) indices, PCN and PCS. PC index values are derived from polar magnetic variations calibrated on a statistical basis such that the index approximate values in units of mV/m of the interplanetary “geo-effective” (or “merging”) electric field (E_M) conveyed by the solar wind. The timing and amplitude relations of the PC index to solar wind plasma and magnetic field parameters are reported. The solar wind effects are parameterized in terms of the geo-effective electric field (E_M) and the dynamical pressure (P_DYN). The PC index has a delayed and damped response to E_M variations and display saturation-like effects for E_M values exceeding 10 mV/m. Steady or slowly varying levels of solar wind dynamical pressure have little or no impact on the PC index above the effects related to E_M for which the solar wind velocity is also a factor. Sharp increases in the dynamical pressure generate impulsive variations in the PC index comprising a initial negative impulse of 5–10 min duration followed by a positive impulse lasting 10–20 min. Typical amplitudes of both the negative and the positive impulses are 0.2–0.5 units. A sharp decrease in the pressure produces the inverse sequence of pulses in the PC index. Auroral substorm activity represented by the AL index level has a marked influence on the average PC/E_M level at the transition from very quiet (AL0 nT) to disturbed conditions while more or less disturbed conditions (AL<−100 nT) have no systematic effect on the average PC/E_M values. At distinct substorm events the PC/E_M ratio has a minimum (0.8) in the pre-onset phase at around 20 min before substorm onset. The average ratio gradually increases in the expansion phase to reach a maximum value (1.1) at around 40 min after substorm onset (or 20 min after the largest (negative) peak in AL). At substorm recovery during the next 2 h the PC/E_M ratio decreases. Finally, we report on the application of polar magnetic variations to model the disturbance storm time (Dst) index development during magnetic storms by using the PC index as a source function to quantify the energy input to the ring current representing accumulated storm energy and characterized by the Dst index.     

Multi-Wavelength Synchrotron Radiation XRF Determination of U and Th in Sedimentary Cores from Lake Baikal

A new technique for the determination of uranium and thorium in lacustrine sediments is based on non-destructive synchrotron radiation X-ray fluorescence analysis (SR-XRF) of sediment samples using monochromatic beams of different energies as excitation sources, with the recording and subsequent iterative processing of the X-ray fluorescence spectra. The technique has a multielement capability and enables fast, simultaneous analysis of a few tens of elements. The proposed SR-XRF technique was tested against INAA and ICP-MS methods and showed a number of advantages with a generally good correlation of results by the three methods. Uranium and Th profiles have been measured at a time resolution of 2 kyr in a drill core (BDP-96) from Lake Baikal bottom sediments deposited between 780 and 40 kyr B P. During this time span, peaks in the U and U/Th concentration marked warm climates. Oscillations of U and U/Th in Baikal bottom sediments during the Brunhes chron reflect climate-induced global change in the volume of polar ice, controlled by long-term cyclicity of the Earths orbital parameters. The response of these warm periods is similar to that observed earlier in shorter cores that span the two last interglacials (220-0 kyr BP).     

Numerical study on transient harbor oscillations induced by successive solitary waves

Tsunamis are traveling waves which are characterized by long wavelengths and large amplitudes close to the shore. Due to the transformation of tsunamis, undular bores have been frequently observed in the coastal zone and can be viewed as a sequence of solitary waves with different wave heights and different separation distances among them. In this article, transient harbor oscillations induced by incident successive solitary waves are first investigated. The transient oscillations are simulated by a fully nonlinear Boussinesq model, FUNWAVE-TVD. The incident successive solitary waves include double solitary waves and triple solitary waves. This paper mainly focuses on the effects of different waveform parameters of the incident successive solitary waves on the relative wave energy distribution inside the harbor. These wave parameters include the incident wave height, the relative separation distance between adjacent crests, and the number of elementary solitary waves in the incident wave train. The relative separation distance between adjacent crests is defined as the ratio of the distance between adjacent crests in the incident wave train to the effective wavelength of the single solitary wave. Maximum oscillations inside the harbor excited by various incident waves are also discussed. For comparison, the transient oscillation excited by the single solitary wave is also considered. The harbor used in this paper is assumed to be long and narrow and has constant depth; the free surface movement inside the harbor is essentially one-dimensional. This study reveals that, for the given harbor and for the variation ranges of all the waveform parameters of the incident successive solitary waves studied in this paper, the larger incident wave heights and the smaller number of elementary solitary waves in the incident tsunami lead to a more uniform relative wave energy distribution inside the harbor. For the successive solitary waves, the larger relative separation distance between adjacent crests can cause more obvious fluctuations of the relative wave energy distribution over different resonant modes. When the wave height of the elementary solitary wave in the successive solitary waves equals to that of the single solitary wave and the relative separation distance between adjacent crests is equal to or greater than 0.6, the maximum oscillation inside the harbor induced by the successive solitary waves is almost identical to that excited by the single solitary wave.     

东昆仑南缘早三叠世洪水川组的源区特征:来自碎屑组成、重矿物和岩石地球化学的证据

分布于东昆仑南缘的早三叠世洪水川组系一套由砾岩、砂岩、页岩和鲕粒灰岩等浅海相和河流相沉积物构成的弧前盆地沉积组合.砂岩碎屑组成、重矿物组合、岩石地球化学和古水流综合研究结果表明,洪水川组砂岩物源主要来自于北侧的岛弧带,同时南侧的阿尼玛卿蛇绿混杂带也为其形成提供了部分物源;源区主要出露长英质岩石、变质岩和硅质岩.     

Water exchange of the Stockholm archipelago—a cascade framework modelling approach

The Stockholm archipelago spans roughly a semicircular area with a radius of approximately 60 km, traditionally partitioned into three parts: the inner, the middle and the outer archipelago. This subdivision coincides with differing water exchange regimes. The inner and middle archipelagos are characterised by comparatively larger basins which are interconnected by a limited number of straits. This configuration is well suited for a discrete basin (DB-) model approach by partitioning the area into a set of sub-basins that are only resolved vertically. The advantage of this approach over 3D-models is the possibility for enhanced vertical resolution and improved strait exchange formulation, outweighing the disadvantage of neglected horizontal gradients within the basins. In the inner archipelago the dominating exchange process is estuarine circulation, induced by the marked freshwater discharge and the vertical mixing. In the outer and middle archipelagos the density fluctuations due to Ekman pumping along the Baltic boundary interface produce another type of baroclinic process that clearly dominates. Measurements to adequately resolve these density variations do not exist. Missing forcing data are provided by linking the middle archipelago's boundary straits to a 3D-model of the Baltic with a grid resolution of 0.5 nautical miles (n.m.). This fine resolution model (FR-domain) is in turn driven by the atmospheric forcing and the density variation at the rectangular boundary of the FR-domain which acceptably resolves both the interfacial straits and the outer archipelago's complex hypsography. Massive computing resources would be demanded if the FR-domain were extended to comprise the entire Baltic. The FR-domain is thus interfaced with an existing coarse resolution model of the entire Baltic (CR-domain) with a grid size of 5 n.m., the open boundary of which is located in the Kattegat. This 3-fold model set-up has been run for one whole year (1992) with a one-year spin-up time to make up for the lack of initial data. The model concept is at this stage to be regarded as a framework for further development in anticipation of improved formulations, particularly for the strait exchange formulation. Therefore only primary validation experiments and a few sensitivity analyses have been performed.